Changes of Interatomic Force Constants Caused by Quantum 
Confinement Eff^ects: Study on the Calculations for the 
First-order Raman Spectrum of Si Nanocrystals in Comparison 

with Experiments. 

Wei-Shan Lee^-^'Q 

^Department of Physics, National Taiwan University, 
Taipei 115, Taiwan, Republic of China. 
^Nano Science and Technology Program, 
Taiwan International Graduate Program, 
Academia Sinica, Taipei 115, Taiwan, Republic of China. 

Abstract 

The redshifts and asymmetric broadening observed in nanocrystal Raman Spectra are attributed 
to the quantum confinement efi'ects by some authors. But others show that they may come from 
the local heating caused by the incident laser as well. In this study we demonstrate that in the 
Si nanocrystal case the latter at most has obvious efi'ects on the broadening but has negligible 
efiFects on the ILO peak shift, while the former contributes most of the ILO peak shift. We 
also demonstrate that after assigning appropriate interatomic force constants in the calculation of 
Raman Spectrum by bond polarizability approximation model within the regime of free boundary 
condition, we may acquire the matching ILO peak shift with experiments. 
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I. INTRODUCTION 



Attributed to its compatibility to the state-of-the-art semiconductor manufacturing pro- 
cess and the abundance of sihcon in nature, sihcon nanocrystals, or in a more favorable 
name, quantum dots, raise world-wide interests for decades in studying their quantum con- 
finement effects on electronic and optical properties'"^ because of the promising applications 
to devices such as the light emitting diodes^ and single electron transistors^. In addition, 
the vibrational properties of the Si nanocrystals also attract great attention^^— . For one 
thing, Raman spectroscopy is a noninvasive technique to investigating defects or qualities of 
devices. For example, the first order Raman shifts may be used to examine the crystalline 
quality of the silicon nanostructures^^^. On the other hand, studies on electron-phonon 
coupling and thermodynamic properties of silicon nanocrystals require full understandings 
for the behaviors of phonon confinements. In comparisons with the first order Raman shifts, 
the second order Raman shifts may be enhanced by electrochemically etched silicon sub- 
strates that are immersed in different duration of tim^. Likewise, Si nanocrystals made by 
annealing hydrogenated amorphous silicon(Q;-Si:II) with continuous-wave laser are also used 
to study the first and the second Raman spectrar^. 

In study of the first-order Raman Spectrum of Silicon nanocrystals, there are at least two 
prominent phenomena to which we should pay attention. First, the asymmetric broadening 
towards low frequency around ILO transition (520 cm"'for bulk) is observed in the spectrum, 
indicating a phonon confinement effect. Secondly, the (bulk) 520 cm"' peak shift to lower 
frequency with reducing nanocrystal sizes. There have been extensive theoretical and exper- 
imental studies of the first-order nanocrystal Raman scattering. H.Richter et.al.— , describes 
a phenomenological exponential function that restricts the nanocrystal phonon wavefunc- 
tions in the sphere. Several different weighting functions are also usedM'—. However, even 
though it is claimed that experimental data could be well explained by the above method, it 
is still obscure for the exponential and weighting functions to have definitive physical mean- 
ings so that they may not provide information of interatomic physical quantities. Therefore, 
it could be more preferable that we start from solving the characteristic vibrations of the 
nanocrystal, which contribute to the change of polarizability resulting in the Raman ef- 
fects. Thus, the bond polarizability approximation!2^(BPA) is introduced to calculate the 
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first-order Raman Spectrum of Si^^ as well as fullereno^. This model suggests that the total 
polarizability tensor is the sum over each axially symmetric polarizability tensor connecting 
each pair of atoms ij, for which may be expressed as a function of two bond- length dependent 
parameters, one referring to mean polarizability while the other anisotropic polarizability. 
The total polarizability tensor is then approximated up to the first order with respect to 
the equilibrium position of the atoms. However, the asymmetric broadening effect does not 
meet agreement with the experimental results in the calculation. 

Failure to success may be imputed to several reasons. First, it is difficult to obtain appro- 
priate interatomic force constants (IFCs) among atoms by only considering the first-neighbor 
interactions, resulting in unsuitable phonon eigenfrequencies and eigenvectors. Improve- 
ments may be made by considering the IFCs up to the second neighbors. Secondly, there 
seems no reason to only consider the first-neighbor atoms in the parameters of bond polariz- 
ability. Thirdly, some studies show that the peak shift and asymmetric broadening may be 
due to local heating from the tiny focusing laser spot rather than size effects^*^. However, 
since both of the parameters of bond polarizability and the IFCs cannot be known in priori, 
it is difficult to obtain proper values of them without comparing the peak shift of exper- 
iments. Therefore, close agreement of asymmetric broadening in concomitant of the peak 
shift towards low frequencies between calculations and experiments remains challenging. 

In this study, we first obtain IFCs by fitting the bulk silicon dispersion relations calculated 
by the rigid-ion model^^ up to the second nearest neighbors with the experimental data from 
the neutron scattering at the three special points, F, X, and L, in the First Brillouin zone. 
The IFCs are then compared with the results from ab initio calculations in order to show the 
significance of the parameters. Afterwards, these fitting parameters are used to calculate 
the vibrational properties in nanocrystals by little modification with the self energy of the 
surface atoms. Nanocrystals with two different diameters are made and measured by Raman 
spectrometer with He-Ne and Argon-ion lasers. Peak shift and broadening effect regarding 
to the nanocrystal sizes or laser power variations are discussed. The first order Raman 
spectra are calculated by bond polarizability approximation (BPA)^. It is found that we 
may acquire close peak shift by assigning proper IFCs in the calculation and obtain more 
satisfactory asymmetric broadening effects by considering the parameters of polarizability 
up to the second-nearest neighbor atoms. Fair resemblance to the calculated first-order 
Raman spectra with the experimental ones is seen. 
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II. THEORY AND SAMPLE PREPARATIONS 



A. The Rigid-ion Model and Interatomic Force Constants in Bulk 



For a homopolar Silicon bulk crystal, the potential energy U is expanded up to the second 
order harmonic approximation with respect to the equilibrium position of an atom. This 



term 



8^U 



^, is then described as the dynamic matrix Djj, in the eigenvalue problem, 



(1) 



a J' 



, where a denotes different kinds of atoms in a unit cell and / denotes x,y,z in Cartesian 
coordinates, being the mass of the atom for each kind. Furthermore, for periodicity in 
bulk, 



(2) 



R 



, where R denotes the position vectors of the first-neighbour or second-neighbour atom and 
S'o-denotes the vectors in which the atom in the primitive unit cell with species a resides. 
In the rigid-ion model^^, the term Djj, {R+ S^^ — 3^^/) is described by several 3 by 3 force- 
constant matrices. Specifically, the 3 by 3 matrix for the first-nearest neighbour atom at the 
position |(1, 1, 1) is defined by 
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and for the second-nearest neighbour atom at the position |(2,2,0) we define 
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. For all the other atoms at the position R, the relative force-constant matrices may be 
derived by similarity transformation (i?) = TiD'^'^ with orthogonal matrix Tj of 

determinant 1^^. Here only up to second-order nearest neighbours are taken into considera- 
tion. We derive the eigenvalues at the F, X, L points in the Brillouin zone^^. For g = F, the 
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eigenvalues are and —8^, both being triplet degenerate. For q = X at (1,0,0), we have 



-AA- 16C ; one LO and one LA (5a) 

-AA-AB-8C-8F ; two TO (5b) 

- + 45 - 8C - 8F ; two TA (5c) 

They are all double degenerate. At last, for the L point (|, |, |), we have 

-6A-2B-8C + 4:D-4:F ; two TO (6a) 

-6A + 4:B-8C -8D-4:F ; one LO (6b) 

- 2A - 4:B - 8C - 8D - AF ; one LA (6c) 

-2A + 2B-8C + 4:D-4:F ; two TA (6d) 



A, B,C, F are solved unambiguously at F and X. Unfortunately one may immediately 
notice that there is no value for D that may simultaneously satisfy Eq 6. Furthermore, 
there is no information for E at the three special points. However, one may obtain values 
of D and E by calculating the interatomic force constants with ab initio calculations. We 
use linear response^*^ implantation of DFT with a local density approximation (LDA) to 
exchange correlation effect and TrouUier- Mart ins pseudopotential in Abinit code^»^. The 
cutoff energy is set to 30 hatrees for the plane-wave basis. Two sets of the parameters 
obtained either by directly solving the multiple equations above, from plugging in neutron 
scattering experiment datar^*^, or by the ab initio calculations are summarized and com- 
pared in Table I. The phonon density of states with even-grid k points^ in the first Brillouin 
zone and phonon dispersion relations are shown in Fig [1] 

One researcher said that the parameter E in the IFCs is zero because "for Group IV, we 
must have i?i = £"2 = to preserve symmetry." While in compounds such as AlAs or GaAs, 
as well as for superlattice, Ei and E2 are also assigned to be zero in their studies. The 
reason why this is incorrect and how they made this mistake will be discussed in Appendix. 

B. Calculations on Nanocrystals Eigenvalue Problems 

An atom in a nanocrystal is called the bulk atom when there are 4 first-nearest neighbours 
and 12 second- nearest neighbours surrounding it, while an atom that does not fulfill this 
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criterion is called the surface atom. Accurate and satisfactory interatomic force constants 
of every atom in a nanocrystal may be obtained by fully considering the minimum value of 
the forces of constraint with arbitrary displacements of atoms in the nanocrystal^^, known 
as the "relaxation method." Nevertheless, with an optimistic simplicity one may place the 
assumption that the interatomic force constants of Silicon atoms in a nanocrystal may 
have close values to those in bulk, except a slight modification at the surface atoms of the 
nanocrystal. The assumption is reasonable because in bulk we only consider that a Si atom 
is influenced at most by its first and second nearest neighbours, leaving outer neighbours 
invisible. Therefore as long as the size of the nanocrystal under consideration is not too 
small, the force constants from bulk may still apply to the nanocrystal. 

At the absence of periodicity, the dynamic matrix of a nanocrystal consisting of atoms 
is expanded to a 3A^ x 3A^ square matrix. The dynamic matrix is constructed in such a way 
that each 3 by 3 block representing the mutual force constants between atom i and atom j 
is placed in positions designating the rows {i — + and columns {j — + 1) in the 
dynamic matrix. The 3x3 diagonal block placed at the rows {i — + 1) and columns 
{i — + 1) in the dynamic matrix is understood as the self energy of the atom i, which 
is discussed below. For bulk, the self-energy of an atom is determined by noticing that the 

net force exerted on an atom at the position {R+ S^) , 

J2 D""' {R + S^-R' - 3^0 ■ u''' (r!) (7) 

R',a' 

, vanishes when the atom is at the equilibrium position d. Thus, 

+ =0 (8) 

R,cr' 

. This gives 

0) = - j2 (R+s.- 3,0 (9) 

. In other words, the self energy of an atom is the negative sign of the sum of the 3 by 3 
force-constant matrices of the atoms surrounding it. In bulk, the self energy is thus 

D'^-(O) = -A{A^ + 2a + F^)I (10) 
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Eq |9] also holds for nanocrystals by modification of the equilibrium position for atoms 
after relaxation. Therefore, for the self energy of surface atoms that have no 16 ambient 
atoms in the nanocrystal we simply sum up the 3 by 3 force-constant matrices of the ambient 
atoms. It may be recognized as the free boundary condition after we make the self energy of 
surface atoms in this way. The central force assumption up to the second-neighbor atoms, 
E = 0^^, is assumed for surface atoms in order to preserve symmetry in the dynamic matrix. 
Afterwards, the standard procedure of solving the eigenfrequencies and eigenvectors of the 
dynamic matrix is applied. Knowledge of group theory regarding to Td group is implemented 
to the dynamic matrix. This dramatically ameliorates the requirements of computer memory 
and reduces time needed to solve the eigenvalue problem of the dynamic matrix. Very few 
number of negative eigenvalues are found because we did not relax the nanocrystals. For 
example, 2 out of 3317 modes in the T2 representation of the 7nm case are found to be 
negative with the imaginary parts 0.494 THz and 0.008 THz. We neglect these kinds of 
modes^. 



C. Bond Polarizability Approximation and Intensity of the Raman Spectrum 

The issue that lies at the heart of the vibrational Raman effects is the change of polariz- 
ability of the nanocrystal on which we measure. The total polarizability tensor where p 
and a refer to the Cartesian coordinate index, is readily understood by summing over every 
polarizability tensor apa{pij) connecting atom i and atom j: 

<* = (11) 

i<j 

Bond polarizability approximation (BPA) model — suggests that one may write the po- 
larizability tensor for each bond apa{pij) as 

(Xpaipij) = a{pij)I + l{pij){pijpij - (12) 

, where a{pij) and '^{pij) are isotropic and anisotropic polarizabilities (called the polariz- 
ability parameters) as a function of position vector associating with atom i and atom j, pij, 
respectively. / is 3 x 3 unit matrix and plj refers to the unit vector of pij. As the nanocrystal 
is excited by incident light the characteristic motion of each atom vibrates with a certain 
eigenfrequency wi away from its own original equilibrium position, therefore the position 
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vector between atom i and atom j changes to pij = r^j + Zjjwith a slight relative move- 
ment Xij ( more specifically, the relative eigenvectors for atom i and atom j corresponding 
to the eigenfrequency wi) displacing away from the relative equilibrium position rjj. For 
harmonic vibrations we may expand Eq [12] with respect to Xij up to the first order under 
the assumption that Xij <^ rij. First notice that, up to the first order. 



A A {Tij + Xij){rij + Xij) 
P^Jp^j = — 



^ij ~r Xij 



~ ''"ij 



^ _ 2(rtj ■ Xij) 



iXij ~^ ■^ij)ij'ij ~^ ■^ij) 



(13) 



Secondly we have the expansion for Oi{pij) up to the first order as 



a{rij + Xij) = a{rij) + a {rij){xij ■ rlj), where a'(rii) ^ 



dxij 



. Similar formula applies to ^{pij)- This leads to the first-order expansion of Eq 
following: 



(14) 

as the 



Oipairij + Xij) ^ a{rij) + a {rij){xij ■ rlj) + -f{rij) + -f' {rij){xij ■ rlj) 



A 



r, ,iy , , _l_ ty , , _l_ rtr> . . ( />■ , 'V • ■ \ \ ■ ■ ■ ■ —I— nf ■ ■ —I— T' ■ ■ I 



a{rij) + a\r,j){xij ■ r%) + 7(ni)(44- " ^) + Vl^jOla^ij ■ 4-)(4-4- - ^) 



Tij ^ ^ Tij ^ , ^ A ^ A A 

" _l_ /y» " I ^ • '/^ 1 "J^ "J^ 



a {rij)I + {rij){rijrij - -) 



^A r)/~^ A,AA 



The first term apcr(?'ij) is neglected because it makes no contribution to Raman effects. Thus 
the total polarizability tensor taken into consideration up to the first-order expansion will 
be expressed as 



a 



tot{ 
pa 



a{nj)I + l\n]){r%r% - 



+7(4)^*/ 



A-^ ^A n,/ A,AA' 

(15) 
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Similarly the expansion of Eq [12] up to the second order may also be derived as 



^(2) ^ 1^3^. AS 
"per 2 ^ ■''^ 



A A 



a {rij)I + {rij){rijrij - -) 



A , A A 



X'ljXij ~h 



A.\ T ■ ■ • T ■ ■ 1 — T . . . — / \ T ■ ■ • V ■ ■ \ \ T ■ T ■ ■ —I— T ■ - T ■ ■ \ 



(16) 



, which involves the second derivative of a and 7 with respect to Xij. In this study, Eq 
will be used to calculate the first-order Raman Spectra for Silicon nanocrystals. 

For backscattering detection configuration, at which the propagation direction of the 
scattered light is parallel but opposite against that of the (linearly polarized) incident light 
, one may obtain the temperature-independent (reduced) intensity of the scattered light 

regarding to two different polarizations: one for the polarization of the scattered light parallel 
to that of the incident light the other for the polarization of the scattered light 

perpendicular to that of the incident light (/s,±)- Taken into consideration of averaging over 
all the possible orientations of nanocrystals for a given radius, we have^^^^: 



(17) 



and 



where giwi) is the vibrational spectrum and 



(18) 



(19) 



a 



tot(l) 
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a 



tot (1)1 2 + [^tot(l) 
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a 



tot (1)1 2 + [^tot(l) 



33 



a 



tot (l)j2" 



tot (1)1 2" 
31 J , 



(20) 



D. Samples Preparations 



The silicon nanocrystals are produced under pressures at 0.5 torr, 1 torr, 2 torr, and 4 
torr, in which we thermal evaporate the silicon sputtering target put on a glass substrate 
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that is mounted on the Tantalum (Ta) boat. The system first is pumped down to 10"^ torr, 
then in tandem, purged with Argon gas three times to remove residual water vapor and 
oxygen. Details of the fabrication process are described elsewhere^^— . In this situation 
the silicon atoms form nanocrystals with various sizes, within some standard deviation, in 
diameters of 7 nm, 13 nm, 31 nm, and 37 nm, respectively. The relationship between the 
fabrication pressure and the size of Si nanocrystals is discussed by elsewhere^^. The standard 
deviations for 7nm and 13nm cases are both roughly Inm. There is a clear tendency that 
under higher pressure leads to a larger nanocrystal size. This may be explained by noticing 
that under higher pressure the higher collision rate is achieved, thus enhancing the migration 
and nucleation probability of the Si nanocrystals in the inert-gas atmosphere. The crystalline 
property is investigated by transmission electron diffraction pattern, and nanocrystal sizes 
are measured by transmission electron microscopy^. 

III. RESULTS AND DISCUSSIONS 
A. Experimental Details 

Raman Spectra of Silicon nanocrystals with two different radii (7 nm and 13 nm in 
diameter) are measured by 488 nm Argon ion laser on Horiba Jobin-Yvon HR800 UV micro- 
Raman Spectrometer with the groove density of 1200 grooves/mm. The charge-coupled 
device detector is cooled down to — 133°C by liquid nitrogen. A lOOA objective is used to 
collect the back-scattering signals. Exposure time of 100 seconds and accumulation number 
of 5 times are required to obtain a clear spectrum in the parallel-polarized configuration for 
every sample. The laser power is varied from 5mW to 60 mW from the source, resulting 
in the laser power on the surface of samples ranging from 0.58mW to 5.23mW, which is 
measured by Thorlabs PMIOOD optical power and energy meter. 

In order to demonstrate the difference in the spectra for polarization dependence, we 
also study the parallel and perpendicular polarizations by 632.8 nm He-Ne laser. Exposure 
time and accumulation number are set to 30 seconds and 20 times and the laser power on 
the sample surface is 1.7 mW in the case. With a resolution of 3.5 cm"^, intensity of the 
polarization-dependent scattered light is measured and compared with the calculation from 
Eq [T7| and Eq [181 A standard Silicon chip, treated as the result of bulk Silicon, is also 
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measured and used to calibrate the system to 520 cm"^. 
B. Calculations of the Spectra 

Before comparing the calculations and experiments, there are some remarks of calcula- 
tions. First, we observe that only intensities involving the eigenmodes belonging to T2, E, 

and Ai are Raman active among the five irreducible representations Ai, A2, E, Ti, and 
T-^. Second, the overall spectrum for a nanocrystal with a specific diameter is acquired 
by summing over intensities from every representation and also by taking into consideration 
the degeneracy of each representation: 

Itot = /ai X 1 + Je X 2 + /r, X 3 (21) 

. Fig EI^A) shows how the three Raman-active representations sum up to a total intensity 
for the case of parallel polarization for the 7nm case. In spite of this, to be consistent with 
the second-neighbor consideration for the short-range part in the dynamic matrix when cal- 
culating the interatomic force constants, we may also need the second-neighbor parameters 
of polarizabilities parameters, a{rij) and 7(rjj). Fig[2](B) shows that the spectrum may be 
matched well for the asymmetric broadening effect for the parallel polarization when one 
considers the polarizability parameters, tt(rjj) and 7(rjj), up to the second-neighbor atoms. 
The parameters of the second-neighbor atoms are chosen to satisfy the ordinary differential 
equations of a(rjj) and 7(rjj)^. This condition, as well as ftc = 7c = 0.0001, is used to 
calculate the first-order Raman spectra of nanocrystals composed of 8,597 atoms (7 nm in 
diameter) and 58,125 atoms (13 nm in diameter). Some anomalous peaks in the range 200 
cm"^ to 300 cm"^ are shown in Fig[2]( 221 cm"^, 254 cm"^, and 281 cm"^ ) but are not shown 
in the experiments because they are mainly due to single-size effects thus being smeared out 
by size deviation in real samples. At last, taking into account the nature of shape of the 
nanocrystals and the fact that the surface of every nanocrystal is hydrogenated, both re- 
sulting in the smaller values of eigenfunctions that stand for the motion of the atoms on the 
surface, we find that the experiments are well described when we divide the eigenfunctions 
in the Ai and E representations, whose eigenvalues are smaller than 300 cm~^, by 8 for the 
7nm case while dividing by 6 for the 13 nm case. 
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C. Comparisons 

1. The Spectra measured by He-Ne Laser 

Comparisons with experiments measured by He-Ne laser with satisfying similarities in 
between are shown in Fig 12] and Fig H] for the polarization-dependent study. The three 
anomalous peaks shown in Fig [2] are smeared out after considering the standard deviation 
(~lnm) of the 7nm nanocrystal diameters. Notice that when using bulk IFCs, the DOS 
of the nanocrystal may almost reproduce that of bulk when the nanocrystal size is larger 
than 6nm, but one immediately recognizes that the mode at 320 cm"^ shown in the DOS of 
nanocrystal does not contribute to the first-order Raman spectrum and that the one-phonon 
effect makes no significant contribution to the lines around 300 cm"^. Therefore it should be 
totally due to the two-phonon effect, being recognized as 2TA— For both cases, the DOS 
curves have the peak at 488 cm"^ while the calculated Raman spectra show the ILO peak 
at 521 cm"^. Because modes in the Ai representation lead to essentially zero intensity in 
Eq [18] , we see that in Fig [3] (B) and Fig H] (B) the spectra with perpendicular polarization 
are flattened out in the 200 cm"^ to 350 cm"^ region. Besides, the peak positions may be 
obtained in calculation by assigning appropriate IFCs. For example, in the 7nm case, we 
may put A=— 9.56431, while in the 13nm case, we may set A = —9.58383. For the rest of 
the IFCs we assign the values with the same ratio of change between bulk value A and the 
assigned new value A. 

To demonstrate the redshifts for smaller nanocrystals, we use the bulk IFCs and compare 
the Raman spectra in Fig O for both polarization directions for nanocrystals with diameters 
1.56 nm and 10.75 nm along with the bulk which shows the delta-like spectrum with peak 
at 524 cm"^. In the inset a clear redshift is shown from 524 cm"^ down to 509 cm"^. One 
observes that the calculated Raman shifts are overestimated ( in the value of 4-7 cm"^) 
compared with those from experiments^^. There are some aspects in this regard. 

The first account for this inconsistency may be that the measured nanocrystals are not 
uniform in size and we should consider along with the standard deviation of the sizes. But 
study shows that the standard deviation for 7nm and 13nm cases is roughly Inm, and in the 
region the Raman shift barely changes according to the inset in Fig [5] This indicates that 
even if we consider the size effects of samples due to standard deviation from the production 
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process it will not improve the peak shifts to get closer to the experimental results. Another 
possible reason may be due to the local heating caused by the laser spot which is roughly 2 
or 3 yum in diameter and has intensity of several milliWatts on the samples^^. Observable 
peak shifts are also claimed to be found when tuning the laser to higher power on the same 
sample. Therefore some studies point out that we may overemphasize the effects of quantum 
confinements on the peak shift. 

2. The Spectra measured by 488nm Argon-ion Laser 

In order to investigate this, spectra changes of parallel-polarized scattered light under 
various incident laser powers in the 7nm and 13nm cases are measured by 488nm Argon-ion 
laser. To obtain the standard deviation, under each laser power the experiment is repeated 
5 times. However, we see that, in Fig|6]^B) and Fig[71^B), the ILO peak shift barely changes 
under different laser powers. This is contradictory to the studies^^^^. The study-i^ points 
out that peak shift influenced by the quantum confinement effect should be negligible for 
nanocrystals larger than 6 nm, which may be shown by our calculation in the inset of Fig 
|5] when using bulk IFCs. However, extrapolation in Fig [6t^B) and Fig El^B) shows that 
under zero laser power the Raman shift should be both around 512 cm"^ for 7nm and 13nm. 
Therefore we believe that the quantum confinement still plays the most important role in 
the redshift (because relaxing the nanocrystal may lead to correct IFCs) for the nanocrystal 
at least with the size up to 13nm. 

Raman peak shifts measured by He-Ne laser for 7nm and 13nm are 513 cm"^ and 516 cm"^, 
respectively. Compared with the zero-power peak shift, 512 cm"^ for both 7nm and 13nm 
cases, they are both within the system measuring error. This indicates that the spectra 
in Fig [3] and Fig S] by 1.7mW He-Ne laser could be thoroughly ascribed to the quantum 
confinement effect; the laser local heating may not contribute in the case. In addition, we 
may infer that the laser local heating may not be very important for the sizes as long as the 
power is lower than 1.7mW. 
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D. Surface Effects on Some Specific Modes 



We now analyze the effects from surface atoms on the three anomalous peaks for the 7nm 
case observed in Fig |5] but not shown in experiments. For comparisons we also include the 
effect on the bulk-like peak, 521 cm"^. To begin with, we study the origin of the four peaks 
by defining the radial probability distribution P\w-wi\<r{i^) for a certain eigenfrequency wi, 
within a natural broadening F, at the distance r away from the center of the nanocrystal 
with the radius R as 



, where Xqj refers to the normalized eigenvector of the atom j. Only modes belonging to Ai, 
E, and T2 are summed in a way similar to Eq[2Il To see the effects from surface atoms we 
may calculate the ratio of integrated probability between surface atoms and bulk atoms: 



, where S and B refer to surface atoms and bulk atoms, respectively. In the case of the 7nm 
nanocrystals, surface atoms refer to atoms in the range of r = 2.26 to r = 2.55. From Fig 
Owe obtain ^^,=521 = 3%, fi^;=22i = 40%, fi^„,=254 = 22%, and fi^,=28i = 33%. Therefore, 
surface atoms barely have an effect on the bulk-like 521cm"^ peak while they may have 
an effect, at most 40%, on the other peaks from single-size contributions. Our samples in 
reahty may be squeezed ellipsoids and this may account for at least 20% why the three 
modes arising partly from surface atoms may not be shown in experiments. Nonetheless, as 
discussed above, the most important factor is that the size deviation in real samples smears 
out the single-size peaks. 

IV. CONCLUSIONS 

In conclusions, we first obtain the interatomic force constants up to the second-neighbor 
atoms, under the rigid-ion model, by directly solving multiple equations that are eigenvalues 
of the dynamic matrix for bulk Si at the three special points F, X, and L in the Brillouin 
zone. The six interatomic force constants (IFCs) thus obtained are compared with the results 




(22) 



A 



(23) 
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from ab initio calculations. Besides matching well for the two sets of parameters, we may 
include as many IFCs for more neighboring atoms as we desire by implementing the results 
from ab initio calculations. The characteristic vibrations of nanocrystals for free boundary 
conditions are solved by assuming that the IFCs among atoms in bulk are optimistically 
applicable to the atoms in the nanocrystals, with a slight modification of the self energy for 

surface atoms. 

Nanocrystals with sizes of 7nm and 13nm in diameter are produced by thermal evap- 
orations and measured in the backscattering configuration on Raman spectrometer with 
632. 8nm He-Ne laser and 488nm Argon-ion laser. The polarization of scattered hght per- 
pendicular or parallel to that of the incident light is distinguished by using an analyzer in 
front of the confocal hole. Besides, the first-order Raman spectra are calculated by bond 
polarizability approximation model. In order to get more satisfactory comparison in the 
spectrum compared with experiments, the two polarizability parameters are extended up 
to the second-neighbor atoms. Under this circumstance, asymmetric broadening is seen in 
fairly good resemblance between calculations and experiments for both sizes measured by 
the He-Ne laser. Moreover, to investigate the redshifts, various sizes of nanocrystals ranging 
from 1.56 nm to 15.36 nm are calculated. We see that when using the same IFCs as bulk 
in the nanocrystals, the peak shifts are almost the same as the diameter of nanocrystals 
is larger than 4nm; the calculated redshifts overestimate the experimental results at the 
amount of 4-7 cm"^, which looks like influenced by the paser local heating effect. But study 
on laser-power dependence show no significant shifts and therefore it is dubious about the 
effect of local heating on the redshifts in the Raman Spectrum. We believe that the local 
heating at most affects the asymmetric broadening on the spectrum shoulders, for higher 
intensity generates more vibrational modes. However, the peak position should not be in- 
fluenced. Therefore we infer that the quantum confinement effect plays the major role in 
the peak shift. Our study shows that we may acquire close peak shift after assigning proper 
values of IFCs whose change is due to the quantum confinement. For example, in the 7nm 
case the value A in IFCs changes to ^4 = —9.56431, while in the 13nm case, we may set 
A — —9.58383. Extrapolation in the 488nm laser-power dependence shows that both 7nm 
and 13nm cases have Raman shift around 512 cm"^ under zero laser power at which the 
peak shift is totally due to the quantum confinement effect. Compared with the shifts in the 
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He-Ne study, 513 cm~^ for the 7nm case and 516 cm~^ for the 13iim case, and noting that 
each of the two values is within the measuring error of the zero-power shift, 512 cm~^, it is 
clear that in the He-Ne study when the laser power is 1.7mW, local heating effect contributes 
nothing to the redshift. 

Finally, the ratios of the integrated probability between surface atoms to bulk atoms 
show that surface atoms contribute only 3% for the bulk-like mode while they make more 
contributions, from 20% to 40%, to the other modes originating from nanocrystals in single 
size. 
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V. APPENDIX 



Some publications^"- set Ei and E2 to be zero, however, they did not give any physical 
explanation. As a matter of fact, one may obtain from the ab initio calculation that neither 
El nor E2 is zero. Indeed, for some materials Ei and E2 are zero due to symmetry property. 
Rocksalt structure, such as NaCl, is an example, but not for the zincblende structure. The 
problem lies in the wrong arrangement in the short-range part of the dynamic matrix when 
looking for the three dimensional force-constant matrices other than that of the atom at the 
position 1(2,2,0) for which we define in EqHl They chose wrong similarity transformation 
Tj, resulting in incorrect symmetry in the short-range part of the dynamic matrix. Every 
element in the dynamic matrix involving Ei and E2 is wrong. The correct elements involving 
El and E2 in the dynamic matrix- for calculating zincblende phonon band structure are 
listed in the row labeled as "correct" in Table II. Terms conjugating to those in Table II are 
also required to be taken into account. However, in the calculations^"^ their counterparts 
are obviously inconsistent. This is due to incorrect assignment for the 3x3 force-constant 
matrices of the second-nearest atoms. 

This mistake seems indifferent for structures with Ei and E2 to be zero, such as NaCl or 
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MgO. But for materials in the diamond structure, such as Sihcon, the mistake may draw to a 
misleading conclusion that "for Group IV, we must have i?i = £'2 = to preserve symmetry." 
In fact, symmetry property in the diamond structure does not restrict the two terms to be 
zero. Kunc, et. al.-, once stated that Ei = E2 = when only considering the central-force 
assumption, but it has nothing to do with symmetry property consideration. Furthermore, 
there is no reason to take Ei = E2 = for materials in zincblende or superlattice, but while 
setting El and E2 to be nonzero, the lethal mistake will give an unpleasant ^^bump" near 
the zone center when exploiting the parameters^ for GaAs, which is shown in Fig. |9] (A). 

To avoid the bump, one may mistakenly replace Ei and E2 to be zero, with the phonon 
band structure shown in Fig. [Hl^B). But it does not give the correct structure for optical 
branches when comparing with the correct ones with appropriate arrangement in the short- 
range part shown in Fig. [9t^C). 

However, it raises a very serious question on whether or not their original arrangement in 
the short-range part produces correct superlattice phonon band structure. Compared with 
the correct short-range part in the dynamic matrix it is easy to find out that the optical 
branches of the superlattice in their publications are essentially incorrect, because they did 
not get correct optical branches for GaAs after assigning Ei and E2 to be zero. Therefore, 
it is imperative to examine again carefully on all of their superlattice-related publications 
with the rigid-ion model. 
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Figures 




r X r L X W LDOS(a.u) 



FIG. 1: Phonon dispersion relation and density of states calculated by Rigid Ion Model 
with the parameters obtained from solving the multiple equations at the three special 
points in the first Brillouin zone as well as E — 0.0189. 
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FIG. 2: (A) The overall spectrum is the sum of the intensities from the three 
Raman-active representations after considering the degeneracy of every representation. 
The curve is for the case with the parameters in the bond polarizability up to the 2nd 
neighbor. (B) Calculations for parallel polarization of the first-order Raman spectrum 
when considering the first-order and the second-order neighboring atoms in the bond 

polarizability parameters. 
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FIG. 3: Comparisons of the first-order Raman spectra between calculations(cal) and 
He-Ne experiments (exp) under parallel (||, Fig A) and perpendicular (±, Fig B) 
polarizations for nanocrystals of 7 nm in diameter. Single-size peaks from 200 cm"^ to 300 
cm"^ are smeared out after summing up all the spectra of the size deviations ranging from 
6.1 nm to 8.4 nm. (cal, ^). After assigning A=-9. 56431, we may get correct peak 
position. Because Intensity from perpendicular-polarized Ai modes is zero, the curve is 
flattened out in the low frequency region in Fig B. 



22 



200- 




60- 



CO 

c 

CD 



40- 



20- 







1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — |- 



T — I — I — I — I — I — I — r 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



(|)=13nm, II 

exp 

cal, S, bulk IPCs 
cal, X, modified IPCs 



(A) 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I M I I I I I I I M m" 



(|)=13nm, ± 

exp 

cal, s, bulk IPCs 
cal, z, modified IPCs 




-i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — i- 

200 250 300 350 400 450 500^ 550 600 650 

Raman Shift(cnn" ) 

FIG. 4: Comparisons of the first-order Raman spectra between calculations(cal) and 
experiments(exp) under parallel (||, Fig A) and perpendicular (±, Fig B) polarizations for 
nanocrystals of 13 nm in diameter. After assigning A=-9. 58383, we may get correct peak 

position. 
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Raman Shift(cm"^) 

FIG. 5: Calculated Raman Shifts vs. nanocrystal diameters with bulk IFCs. The 1.56 nm 
nanocrystal shows the peak at 509 cm"^ while the 10.75 nm one at 521 cm"^. Bulk has a 
shift of 524 cm"^. (Inset) The calculated results, shown in square, overestimate the 

experimental ones, shown in circle. 
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FIG. 6: Power variations by 488nm Argon-ion laser on the 7nm case. (A) Raman Spectra 
under different laser powers on the sample surface. (B) ILO Raman Shift vs. laser power 
on the sample surface. No significant shift is observed under different laser powers. 
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FIG. 7: Power variations by 488nm Argon-ion laser on the 13nm case. (A) Raman Spectra 
under different laser powers on the sample surface. (B)ILO Raman Shift vs. laser power 
on the sample surface. No significant shift is observed under different laser powers. 
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FIG. 8: Radial probability distributions for the 7nin case. (A) Bulk- like Raman shift. 
Prom (B), (C), and (D) we calculate that the ratios of integrated probability are 40%, 

22%, and 33%, respectively, a = 5.43. 
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FIG. 9: GaAs phonon band structure. (A) When setting nonzero Ei and E2 in their 
original fortran procedure, an unpleasant bump is shown near the zone center owing to the 
wrong symmetry in the short-range part in the dynamic matrix. (B) When assigning the 
El and E2 to be zero, the unwanted bump is removed, but we have incorrect optical 
branches while comparing with the correct band structure in (C). 
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Tables 



TABLE I: IPCs obtained either by solving Eq. 5 and Eq. 6 with neutron scattering data or 



by Abinit calculation. 





A 


B 


C 


F 


D 


E 


Neutron Scattering Data"^ 


15.69 




13.69 
12.45 
4.51 




14.62 
12.51 
10.97 
3.41 




Solutions from Eq. 5 and Eq. 6* 


-9.8315'= 


-6.6725 


-0.6373 1.4045'' 


-0.5210^= 


0.0189^ 


Abinit calculations'' 

'J- „i' '|1U„ 


-9.1264 


-6.4525 


-0.4970 


1.1317 


-0.4781 


0.2944 



^In unit of kg/sec^ after each parameter value is multiplied by 5.7641. 
^At r. 

''At X by solving Eq. 5. 

^At L. D \s chosen to get the closest fit to Eq. 6. 

is chosen to fit the values in the L-X-W-L direction. 



TABLE IL Comparisons between correct and wrong elements involving the parameters Ei 
and E2 in the dynamic matrix for the short-range part. 
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correct 
wrong 


— 4I?i sin(7rA;i^ 
—ADi sin(7rA:i^ 


sin(7rA;2) 
sm(7rA;2) 


+ AiEi sm(7rA;3) cos(7rA;i^ 
+ AiEi sin(7rA;i) cos(7r/c3^ 


— iiEl sin(7rA;3 

— AiEi sin(7rA;2 


cos(7rA;2) 
cos(7rA;3) 




correct 

wrong 


—AD\ sin(7rA;i^ 

— IDi siii(7r/,:i' 


sinfTr/uQ 

siu(7r/,-:-!) 


+ AiEi sin(7rA;2) cos(7rA;i' 

+ l/Ei siu(7r/,-i) cos(7r/,'2) 


— AiEi sin(7rA;2 

— UEi sm(TTk:',] 


cos(7rA;3) 

C<)s(^/,-2) 




correct 

wrong 


—ADi sin(7rfc2; 

—ADi sin(7r/c2^ 


sin(7rfc3" 

sm(7rA;3) 


+ 4iEi sin(7rfci) cos(7rA;2' 

+ AiEi sm(7r/c2) cos(7r/ci^ 


— 4iE\ sin(7rA;i' 

— 4iEi sin(7rA;3 


cos(7rA;3) 

cos(7rA;i) 




correct 
wrong 


— 4D2 sin(7r/cij 
— 4Z)2 sin(7rA;i; 


sin(7rA;2 
sin(7rA;2^ 


— 4i£'2 sin(7rA;3) cos(7rfci' 
+ 4iE2 sin(7rA;i) cos(7rA;3] 


+ 4i£^2 sin(7rfc3 
— AiE2 sin(7rA;2^ 


cos (77^:2) 
cos(7rA;3) 




correct 
wrong 


— 4Z)2 sin(7rA;i; 
— 4D2 sin(7r/ci^ 


sin(7rfc3" 

sin(7r/c3^ 


— 4iE2 sin(7rA;2) cos(7rA;i' 
+ 4i£?2 sin(7rA;i) cos(7r/i:2' 


+ AiE2 sin(7rA;2^ 
— 4iE2 sin(7rfc3' 


cos(7rA;3) 

cos(7rA;2) 




correct 
wrong 


— 4D2 sin(7r/c2^ 
— 4I?2 sin(7rA;2; 


sin(7rA;3) 
sin(7r/e3) 


— 4i£^2 sin(7rA;i) cGs(7rA;2^ 
+ 4i£^2 sin(7r/c2) cos(7rA;i) 


+ 4i£'2 sin(7rA;i 
— 4i£'2 sm.{'Kk2>\ 


cos(7rfe3) 
cos(7rA;l) 
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